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Abstract 

The eigenvectors of the osp(l\2) invariant Gaudin hamiltonians are found 
using explicitly constructed creation operators. Commutation relations between the 
creation operators and the generators of the loop superalgebra are calculated. The 
coordinate representation of the Bethe states is presented. The relation between the 
Bethe vectors and solutions to the Knizhnik-Zamolodchikov equation yields the norm 
of the eigenvectors. 
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1 Introduction 



Classifying integrable systems solvable in the framework of the quantum inverse scattering 
method Q by underlying dynamical symmetry algebras, one could say that the Gaudin 
models are the simplest ones being related to loop algebras and classical r -matrices. 
More sophisticated solvable models correspond to Yangians, quantum affine algebras, 
elliptic quantum groups, etc. 

Gaudin models || are related to classical r -matrices, and the density of Gaudin 
hamiltonians 

H ia) =J2 r ab(Za-Zb) (1-1) 

coincides with the r -matrix. Condition of their commutativity [H^ a \ H^] = is 
nothing else but the classical Yang-Baxter equation (YBE). 

The Gaudin models (GM) associated to classical r -matrices of simple Lie algebras 
were studied in many papers (see [3-9] and references therein). The spectrum and eigen- 
functions were found using different methods (coordinate and algebraic Bethe Ansatz 
II ||, separated variables ||, etc.). A relation to the Knizhnik-Zamolodchikov (KZ) 
equation of conformal filed theory was established [7-9] . 

There are additional peculiarities of Gaudin models related to classical r -matrices 
based on Lie superalgebras due to Z2 -grading of representation spaces and operators. 
The study of the osp{l\2) invariant Gaudin model corresponding to the simplest non- 
trivial super-case of the osp{\\2) invariant r -matrix flCi[| started in |TT|| . The spectrum of 
the osp(l\2) invariant Gaudin hamiltonians was given, antisymmetry property of 

their eigenstates was claimed, and a two site model was connected with some physically 
interesting one (a Dicke model). 

The creation operators used in the sl(2) GM coincide with one of the L -matrix 
entry ||, |4j. However, in the osp(l\2) case, as we will show, the creation operators are 
complicated polynomials of two generators X + (X) and v + (/j,) of the loop superalgebra. 
We introduce B -operators by a recurrence relation (Section 3). Acting on the lowest spin 
vector the B -operators generate exact eigenstates of the Gaudin hamiltonians , 
provided Bethe equations are imposed on parameters of the states. Furthermore, the 
recurrence relation is solved explicitly and the commutation relations between the B- 
operators and the generators of the loop superalgebra C(osp(l\2)) are calculated. We 
prove that the constructed states are lowest spin vectors of the global finite dimensional 
superalgebra osp{\\2) , as it is the case for many invariant quantum integrable models 
[j~2[ ]. Moreover, a striking coincidence between the spectrum of the osp(l\2) invariant 
Gaudin hamiltonians of spin s and the spectrum of the hamiltonians of the sl(2) GM 
of the integer spin 2s is found (Section 3). 

A connection between the B -states, when the Bethe equations are not imposed on 
their parameters, of the Gaudin models for simple Lie algebras to the solutions to the 
Knizhnik-Zamolodchikov equation was established in the papers jj], ||. An explanation 
of this connection based on Wakimoto modules at critical level of the underlying affine 
algebra was given in ||. An explicit form of the Bethe vectors in the coordinate repre- 
sentation was given in both papers jj], ||. The coordinate Bethe Ansatz for the B -states 
of the osp(l\2) GM is obtained in our paper as well. Using commutation relations be- 
tween the B -operators and the transfer matrix t(X) we demonstrate algebraically that 
explicitly constructed B -states yield a solution to the KZ equation (Section 4). This con- 
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nection permits us to calculate the norm of the eigenstates of the Gaudin hamiltonians. 
An analogous connection is expected between quantum osp(l|2) spin system related to 
the graded YBE [10, 13] and quantum KZ equation following the lines of JUp. 



The norm and correlation functions of the sl(2) invariant GM were evaluated in 
H using Gauss factorization of a group element and Riemann-Hilbert problem. The 
study of this problem for the GM based on the osp(l\2) Lie superalgebra is in progress. 
Nevertheless we propose a formula for the scalar products of the Bethe states which is 
analogous to the sl(2) case (Conclusion). 

2 OSp(l\2) -invariant i?-matrix 

Many properties of the Gaudin models can be obtained as a quasi-classical limit of the 
corresponding quantum spin systems related to solutions R(X; rf) to the YBE. In the 
quasi-classical limit ij — > 

R(X;rj) = I + 7 1 r(X) + 0(7 1 2 ), 

some relations simplify and therefore can be solved explicitly providing more detailed 
results for the GM. 

The graded Yang-Baxter equation || differs from the usual YBE by some sign 
factors due to the embedding of R -matrix into the space of matrices acting on the Z2 - 
graded tensor product V\ ® V% <S> V3 . At this point our aim is to reach fundamental 
osp(l\2) invariant solution. The rank of the orthosymplectic Lie algebra osp(l\2) is one 
and its dimension is five. The three even generators are h, X + ,X~ and the two odd 
generators are v + ,v~ |L7]]. The (graded) commutation relations between the generators 
are 

[h,X ± ] = ±2X ± , [X+,X-]=h, 
[h,v ± }=±v ± , [v+,v-] + = -h, (2.1) 
[X*,v ± ] = , [v ± ,v ± } + = ±2X ± , 

together with [A ± ,f ± ] = . Notice that the generators h and considered here fl2.1|) 



differ by a factor of 2 from the ones used in [17, [H]]. Thus the Casimir element is 



c 2 = h 2 + 2 (X + X~ + X~X+) + (v + v- - v~v + ) 

= h 2 - h + 4A+ X~ + 2v + v- . (2.2) 

For further comparison with the sl(2) Gaudin model |3|, [| and due to the chosen set of 
generators ( [2.1J ) we parameterize the finite dimensional irreducible representations V^' 
of the osp(l|2) Lie superalgebra by an integer /, so that their dimensions 2/ + 1 and the 
values of the Casimir element ( |2.2D C2 = /(/ + 1) coincide with the same characteristics 
of the integer spin I irreducible representations of sl(2) . 

The fundamental irreducible representation V of osp(l\2) is three dimensional. We 
choose a gradation of the basis vectors e^; j = 1, 2, 3 to be (0, 1, 0) . 



The invariant R -matrix is a linear combination [1C] 



R = X(X + ^-)l + r ] (x + ^-)V-7 ] XK, (2.3) 



2 'V 2 



of the three OSp(l\2) group invariant operators [g (g> g, X] = , g € OSp(l\2) , X € 
End (V (g> V) , acting on V ® V : the identity / , the graded permutation V and a 
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rank one projector K . In the equation (2.3) A is the spectral parameter, and r\ is a 
quasi-classical parameter. 

The L -operator of the quantum spin system on a one-dimensional lattice with N 
sites coincides with R -matrix acting on a tensor product Vq V a of auxiliary space Vq 
and the space of states at site a = 1,2, ... N 

L 0a {\ - z a ) = R 0a {\ - z a ) , (2.4) 

where z a is a parameter of inhomogeneity (site dependence) [l], |2| . Corresponding mon- 
odromy matrix T is an ordered product of the L -operators 

N 

T(X; {z a }f ) = L 0N (X -z N )... L 01 (X - z x ) =[J L 0a (X - z a ) . (2.5) 

a=l 

The commutation relations of the T -matrix entries follow from the FRT-relation Q 

Ri2{\ - m)Ti(A)T 2 (m) = r 2 (/x)Ti(A)i?i2(A - n) . (2.6) 



Multiplying ( |2.6[) by i2 12 x and taking the super-trace over V\ (g) V2 , one gets commuta- 
tivity of the transfer matrix 

*( A ) = E< l)'- ; V),iA: {za}?) = T U - T 22 + T 33 (2.7) 

for different values of the spectral parameter t(X)t(/j,) = t(fi)t(X) . 

The choice of the L -operators ( p.4| ) corresponds to the following space of states of 
the osp(l\2) -spin system 

N 

n = ® v a . 

a=l 



The eigenvalue of the transfer matrix i(A) in this space is |1C] 

M 

A(A;{/, J }f) = a{ N) (X;{z a }f)Hs l (X- H )-ai N) (X;{z a }f)x 

3=1 

M 

x Jj5i(A-/i i + |)5_i(A-/4 i + r7) + 

+ 4 N) (A;Wf)n^(A- ft + ^, (2.8) 

where aj^ (A; {z a }? ) = n&Li ay (A - z 6 ) ; j = 1,2,3, 

ai(A) = (A + 77) (A -+- 3^/2) , a 2 (A) = A(A + 3? ? /2) , 

a 3 (A) = A (A + r,/2), S n (jM) = ^~^ /2 2 . (2.9) 



Although according to (|2.8|) the eigenvalue has formally two sets of poles at A = nj — rj/2 
and X = /ij — r] , the corresponding residues are zero due to the Bethe equations [|n]] 

N , _ /9 v A/ 

n - ,0 ) = n ^ - ^)s-2in - m*) • (2.10) 
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If we take different spins l a at different sites of the lattice and the following space of 
states 

N a \ 
H=® V} la > , 

a=l 

then the factors on the left hand side of (|2.10| ) will be spin dependent too 

fJ-j ~ z a + rjlg/2 
fJ-j ~ z a - r]l a /2 ' 

The osp{\\2) invariant R -matrix (|2.3|) has more complicated structure than the sl{2) 
invariant R -matrix of C. N. Yang R = XI + rfP . As a consequence the commutation 
relations of the entries Tij(X) of the T -matrix (|2.5|) are more complicated and construc- 
tion of the eigenstates of the transfer matrix t(X) by the algebraic Bethe Ansatz can be 
done only using a complicated recurrence relation expressed in terms of Tij(/j,k) (l6| (see 
also for the case of osp(l\2) ). It will be shown below that due to a simplification 
of this relation in the quasi-classical limit ij — > one can solve it and find the creation 
operators for the osp{\\2) Gaudin model explicitly. 



3 OSp(l 1 2) Gaudin model and creation operators 



The classical r -matrix of the orthosymplectic Lie superalgebra osp(l\2) can be expressed 
in a pure algebraic form using Casimir element in the tensor product osp(l\2) (g> osp{\\2) 





f(A) 



-2 ) 



(3.1) 



where cf = h®h + 2 (X + ® X~ + X~ X + ) + (v + ® v~ - v~ <8> v + ) . The matrix 
form of the Casimir element f in the fundamental representation n of osp{\\2) follows 
from ( |3.l[) by substituting appropriate 3x3 matrices instead of the osp(l\2) generators 
(2.1) and taking into account ^-graded tensor product of even and odd matrices [10|. 
Alternatively, the same matrix form of r can be obtained as a term linear in ij in the 
quasi-classical expansion of ( glj ) 

r(\) = r j-= l -(V-K) , 

where V is a graded permutation matrix and K is a rank one projector. 

A quasi-classical limit r\ — > of the FRT-relations ( |2.6| ) results in a matrix form of 
the loop superalgebra relation (T(A; rj) = I + r]L(X) + C(r/ 2 )) 



L (A), L (m) 



ri 2 (A-p), L (A)+ L (/i) 



(3.2) 



Both sides of this relation have the usual commutators of even 9x9 matrices L (A) = 

L(A) <S> I3 , L (a*) = ^3 ^(m) and ri 2 (A — /x) , where I3 is 3 x 3 unit matrix and L(A) 
2 

has loop superalgebra valued entries: 



L(A) 



MA) 
«+(A) 
2A+(A) 



-«-(A) 2X-(X) 
u-(A) 
u+(A) -/i(A) 



(3.3) 
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The relation ( |3.2| ) is a compact matrix form of the following commutation relations 
between the generators h(X), v ± (fi), X^(v) of the loop superalgebra under consideration 



[h(X) , X ± (fi)] = T2 
[X+(X) , X-Qjl)] = 
[v+ (A), v-fr)] .= 



X*(A)- **(//) 
A — /i 
_ /t(A) - 
A — 
fe(A) - /t(/x) 
A — /i 



[h(X) , V ±( M )] = 
[X±(A) , ^( M )] 
[^(A), ^(m)] . 



u*(A) - u^Qf) 

X — fj, 
v ± (X) - v ± (fi) 

X — fj, 
X±(A)-AT±( M ) 



A — /i 



(3.4) 



together with [X ± (A) , w ± ( / u)] =0. 

These commutation relations (3.4) define the positive part £ + (osp(l\2)) of the loop 
superalgebra. The usual generators Y n of a loop algebra parameterized by non-negative 
integer, are obtained from the expansion 

Y ( x ) = E„>o Y n A-(" +1 ) . In particular, taking 
all Y n = for n > one gets an L -operator L(A) = Lq A -1 , where Lo is osp(l\2) - 
valued matrix. This Lq similar to tq satisfies cubic characteristic equation with the 
osp(l|2) Casimir element (|2.2| ) as coefficient 

Ll + 2L\ - (c 2 - 1) L Q - pal = . (3.5) 

A Gaudin realization of the loop algebra (|3.4| ) can be defined through the generators 
y = (/ l ,t; ± ,X ± ) 



N 



YW = E 



a=l 



X- Z a 



Y a G End (K) , 



(3.6) 



where are osp(l|2) generators in an irreducible representation V a of the lowest 
spin —l a associated with each site a || ||. Then the L -operator (|3.3| ) has the form 



^(A;K}f) 



N 

£ 

0=1 



(3.7) 



here {z a }^ are parameters of the model (cf ([2.5|)). It follows from the relation (|3.7| ) that 
the first term in the asymptotic expansion near A = oo defines generators of the global 
superalgebra osp(l\2) C C + (osp(l\2)) 



N 



L 



gl 



lim AL(A) 

A^oo 



(3. 



a=l 



with h g i , v i and Af ^ as generators and entries of L g \ (cf ( |3.3| )). Moreover, from the 
equation (|3.2[) we get 



. l 2 

here L 9i = L g ; <g> I 3 , e.#. [/i^ , w + (/u)] = v + {p) . 



ro , L (m) 
2 



(3.9) 
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Let us consider the loop superalgebra £+(osp(l|2)) as the dynamical symmetry al- 
gebra, i.e. as the algebra of observables. In order to define a dynamical system besides 
the algebra of observables we need to specify a hamiltonian. It is a well-known fact that 
due to the r -matrix relation (|3,2| ), the so-called Sklyanin linear brackets, the elements 



t(X) = ±stTL 2 (\) = h 2 (\) + 2[X+(\),X-(\)} + + lv + (\),v-(\)}- 

= h 2 {\) + ti(\) + 4X + (\)X-{\) + 2v + (\)v-(\) (3.10) 

commute for different values of the spectral parameter t(X)t(fi) = t(fi)t(X) . Thus, t(A) 
can be considered as a generating function of integrals of motion. 

It is straightforward to calculate the commutation relations between the operator 
t(A) and the generators of the loop algebra X + (p>) and v (/x) 

[t(A),X+Qu)] - j X+ M h W ~ X+ W h ^) v + (X)v + (p)-v+{p)v + (X) ^ (3 n) 

A — fi A — (J 

[t(X),v+M] =2 » + WA)-^WW _ 4 X + („)„-(A)-X + (A)„- M 

A — /i A — 



A direct consequence of the equation ( |3.9| ) is an invariance of the generating function 
of integrals of motion t(X) under the action of the global osp(l\2) [t(X),L g i] = . 

We can consider the representation space H p h of the dynamical algebra to be a lowest 
spin ^(A) representation of the loop superalgebra with the lowest spin vector f2_ 

h(X)Q- = p(X)Q- , u~ (A)Q_ = . (3.13) 



In particular, a representation of the Gaudin realization (3^7) can be obtained by con- 
sidering irreducible representations Va of the Lie superalgebra osp{\\2) defined by a 
spin —l a and a lowest spin vector uj a such that v~u a = and h a u a = —l a ^ a ■ Thus, 

N X - —In 

Q_= ® Ua and p(X) = > — . (3.14) 

a=l , X- Z a 

It is a well-known fact in the theory of GM |j, that the Gaudin hamiltonian 

F(a)= y-£|M) ) (315) 

here cf (a, 6) = /i a /i(, + 2 (X+X^ + X^X^) + (v+v^ — v~v^) , can be obtained as the 
residue of the operator t(X) at the point A = z a using the expansion 



0=1 



To construct the set of eigenstates of the generating function of integrals of motion 
t(X) we have to define appropriate creation operators. The creation operators used in 
the si (2) Gaudin model coincide with one of the L -matrix entry |3|, |4j. However, in 
the osp(l\2) case the creation operators are complicated functions of the two generators 
X + (A) and v + (p) of the loop superalgebra. 
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Definition 3.1 Let Bm((J>i, ■ ■ ■ , Mm) belong to the Borel subalgebra of the osp(l\2) loop 
superalgebra £+(osp(l\2)) such that 

M (-IV 

B M (fj, 1 , . . . ,(i M ) = v + (fj, 1 )B M ^ 1 (/i 2 , . . . ,/iAf) + 2X + Oi) y~] Btyjfja,.. .,hm) , 

(3.17) 

with Bq = 1, Bi(fi) = and Bm = for M < . The notation 

B^j_ 2 (/j,2, ■ ■ ■ means that the argument jjLj is omitted. 

As we will show below, the B -operators are such that the Bethe vectors are generated 



by their action on the lowest spin vector f2_ (|3.13[) , To prove this result we will need some 



important properties of the B -operators. All the properties of the creation operators 
Bm(^i, ■ • • j Mm) listed below can be demonstrated by induction method. Since their 
proofs are lengthy and do not contain illuminating insights we will omit them. 

1 The creation operators £?m(Mi> • • • j Mm) are antisymmetric functions of their argu- 
ments 

Bm(pi, ■ ■ ■ , Vk, Mfc+ij • • • j Mm) = - Bm(hi, ■ ■ ■ ,Hk+i,Hk, ■ ■ ■ >Mm) , 

(3.18) 

here 1 < k < M and M > 2 . 

Subsequently we calculate the commutation relations between the generators of the 
loop superalgebra £+{osp{\\2)) and the B -operators. In order to simplify the formulas 
we will omit the arguments and denote the creation operator f?M (Mil ■ • • > Mm) by Bm ■ 

2 The commutation relations between the creation operator Bm and the generators 
v + (X) , h(\) , v~(X) of the loop superalgebra are given by 

v+(X)B M = (-l) M B M v + (X) + 2 X>l) i X+(A ,~ X+(Mj) B M l, , (3.19) 



■ w 1 \ £ (-1)* 



h{X)B M = B M h(X) + ]T y— + 



x 



Mi / ~t^~Mi 

,.v M f_lV'+ e (*-J') 



(^)- B M-i + 2A+(A) £ V /._„, Bj& , (3-20) 



, , y A "Mj ~ (A-Mfe)(Mfe -Mi) 

M D (i,j) 



Mi ~ Mj V A ~ Mr A ~ Mj / 

(3.21) 

(7) 

/iere £/ie upper index of B M-1 means that the argument fj,j is omitted, the upper index 
of bH_ 2 means that the parameters are omitted and is Heaviside function. 
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Already at this point we can make some useful observations. 

Remark 3.1 The commutation relations between the generators of the global osp(l\2) 
and the B -operators follow from the previous property. To see this we multiply ( 3.1$) , 
(3.2C) and (3.21) by A and then take the limit A — > oo . In this way we obtain 

M 

v+B M = (-fBM.j-SEt-^rf-i- ( 3 - 22 ) 
h gl B M = B M (h g i + M) , (3.23) 

M I M 7 

%iBm = (-i) M B M v; l + Y,(- i y B M-A h ^) + E^— l 

(3.24) 

The proof of the main theorem is based on subsequent property of the creation operators. 

3 The generating function of integrals of motion t(X) j[3.1(\ ) has the following commu- 
tation relation with the creation operator Bm{^i, ■ ■ ■ ,}J>m) 



I M 1 - 1 
t(X)B M = B M t{\) + 2B M h(X) V + V -r r- 




£ (X - Ih)(X - fij) 

M (—^ 

+ 4^f^< ) _ 1 (l + (AK( /Jl )-I+( /l! K(A)) . 
i=i A " ^ 

XYie notation we use here for the operator fluim) = h(m) + J^f^ii^i ~ Mj) 1 • 

In the Gaudin realization ( |3.6| ) the creation operators £?m(m1j ■ ■ ■ , ^m) have some 
specific analytical properties. 

4 The B -operators in the Gaudin realization ( \3.6[ ) satisfy an important differential 
identity 

This identity will be fundamental step in establishing a connection between the Bethe 
vectors and the KZ equation. 

The recurrence relation ( 3. 17[) can be solved explicitly. To be able to express the 
solution of in a compact form it is useful to introduce a contraction operator d . 

Definition 3.2 Let d be a contraction operator whose action on an ordered product 

M 

\\ v + (fij ) , M > 2 , is given by 

3=1 

M-l M ( _, )(T (jk) M 

d{v + ( fH )v + ( t i 2 )...v+( m ))=2j2x + (H) E II v+ ^ml (3.27) 

j=l k=j+l ^ k m^j,k 
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where cr(jk) is the parity of the permutation 

a : (l,2,...,j,j + l,...,fc,...,M) -» (1, 2, . . . , j, k, j + 1, . . . , M) . 

M 

The d operator can be applied on an ordered product J\ v + (fj,j) consecutively several 

3=1 

times, up to [M/2] , the integer part of M/2 . 

Theorem 3.1 Explicit solution to the recurrence relation ( 3.1TQ is given by 



M [M/2] M 

. • ■ ,w) =n + e — dm n v+ ^) ■ ( 3 - 28 ) 



3=1 771=1 3=1 



The properties of the creation operators Bm studied in this Section will be funda- 
mental tools in determining characteristics of the osp(l\2) Gaudin model. Our primary 
interest is to obtain the spectrum and the eigenvectors of the generating function of 
integrals of motion t(A) ( |3,10| ). 



Theorem 3.2 The lowest spin vector S7_ (3.13J is an eigenvector of the generating 



function of integrals of motion i(A) ([3. 1 (\ ) with the corresponding eigenvalue Ao(A) 

t(A)fi_ = A (A)O_ , A (A) =p 2 (X)+p'(\) . (3.29) 

Furthermore, the action of the B -operators on the lowest spin vector 0_ yields the 
eigenvectors 

^(m, ■ ■ ■ , hm ) = B M (m, . . .,Hm) ^- , (3.30) 

of the t(X) operator 

*(A)*( Mi ,...,mm)=A(A; {ii j }f =l )^{ii 1 ,... 1 iiM) , (3.31) 
with the eigenvalues 

M 

A(A; { N }f =l ) = A (A) + 2p{\) £ + 2 £ — , (3.32) 

provided that the Bethe equations are imposed on the parameters {^j}jL\ 

M 

fabi) = p(h) + Yl — = • ( 3 - 33 ) 



Proof. The equation (|3.29|) can be checked by a direct substitution of the definitions 
of the operator i(A) and the lowest spin vector 0_ , the equations ( p. 10 ) and ( 3.13|) , 
respectively. 



To show the second part of the theorem, we use the equation ( 3.30| ) to express the 
Bethe vectors ^(pi, • • • , Pm) 

t(\)y(fn,...,n M )=t(\) n- • (3.34) 



9 



Our next step is to use the third property of the B -operators, the equation ( |3.25| ), and 
the definition of the lowest spin vector S7_ ( |3.13| ) in order to calculate the action of the 
operator t(X) on the Bethe vectors when the Bethe equations ( |3. 33 ) are imposed 



/ M 

t(\)B M Sl- = B M t(X)n- + 2 I p(X) Y — 



M 



1 



fM + ^ (A -//,!( A -//,) 



BmO- . (3.35) 



We can express the first term on the right hand side since we know how the operator 
t(A) acts on the vector J7_ , the equation (|3.29|) , 



t(\)B M n_ 



Ao(A) + 2 



(M M 
i=i Pl i<j 



B M ^- . (3.36) 



The eigenvalue equation ( |3.31| ) as well as the expression for the eigenvalues ( |3.32 ) follow 
from the equation ( 3.36|) . □ 



Corollary 3.1 In the Gaudin realization of the loop superalgebra given by the equations 
( \3. (\) and ( 3.14 ) the Bethe vectors ^(/^i, • • • , Hm) ( S.SCj) are the eigenvectors of the 
Gaudin hamiltonians ( 3.1$ ) 



with the eigenvalues 



E 



(a) 
A I 



N 



■y- — y.-L 



M 



la 



b^a 



Zb 



=1 z a 



(3.37) 
(3.38) 



when the Bethe equations are imposed 

M 



pm(m>j) = p(j*j) + Y 



1 



V 



M 



fJ'j - Mfc 



E-^r + E 



1 



a=l 



k+j 



(3.39) 



The statement of the corollary follows from residue of the equation ( |3.3l| ) at the point 
A = z a . The residue can be determined using ( |3,16j ), ( |3.32 ) and ( 3.29[ ). 

The eigenvalue ( p.32| ) of the operator t(X) and the Bethe equations ( |3.33| ) can be ob- 
tained also as the appropriate terms in the quasi-classical limit rj — > of the expressions 
Q and ( plO| ). 



Comparing the eigenvalues ( 3.38] ) of the Gaudin hamiltonians and the Bethe 

equations ( |3,39| ) with the corresponding quantities of the sl{2) GM |Ij we arrive to 
an interesting observation. 

Remark 3.2 The spectrum of the osp(l\2) Gaudin model with the spins l a coincides 
with the spectrum of the sl(2) GM for the integer spins (cf. an analogous observation 
for partition functions of corresponding anisotropic vertex models in [HjJ). 



Remark 3.3 The Bethe vectors are eigenstates of the global generator h g i 



(3.40) 
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Moreover, these Bethe vectors are the lowest spin vectors of the global osp(l\2) since 
they are annihilated by the generator v~ L 

v- l V(n 1 ,..., m ) = 0, (3.41) 

once the Bethe equations are imposed ( 3. 3!\ ). These conclusions follow from the remark 
3.2 the equations ( 3. 2c ), ( 3.2^ ) and the definition of the Bethe vectors ( 3. 3C ). 

Hence, action of the global generator on the lowest spin vectors ^(//i, ••• 
generates a multiplet of eigenvectors of the operator t(A) . One can repeat the arguments 
of HI to demonstrate combinatorially completeness of the constructed states. 

As was pointed out already in || for the sl(2) case, there are several modifications 
of the hamiltonians ( 3. 15| ) . One of them is the Richardson's pairing- force hamiltonian. 
These modifications can be formulated in the framework of universal L -operator and 
r -matrix formalism ( |3.2| ) Q. 

Due to invariance of the r -matrix (34) one can add to the L -operator any element 
of osp(l\2) 



L(A) - L(A) = gY + L(A) , 
preserving commutation relations ([T^) • If we choose Y = h , then 
t(X) = \ str L 2 (A) = i(A) + 2g h(X) + g 2 



(3.42) 



(3.43) 



will have the commutativity property, i.e. t(X)t(fj,) = t(fj,)t(X) . Hence we can take t(X) 
to be the generating function of the (modified) integrals of motion 



cf(a,b) 



Zb 



(3.44) 



Notice that the global osp(l\2) symmetry is now broken down to global u(l) . In this 
case the eigenstates \£ are generated by the same B-operators, but eigenvalues and 
Bethe equations are slightly changed. Richardson like hamiltonian || can be obtained 
as a coefficient of 1/A 2 in the A — > oo expansion 



t(X) 



2 i 2 9 u 

9 + y h g i + A2 



^ (^2g J2 z ah a + c 2 (gl)j +0(^; 



(3.45) 



Another modification can be obtained using an L -operator with bosonic and fermionic 
oscillator entries 

fX - 7 2b \ 
Losc(X) = 7 + 7 , (3.46) 
\ 26+ 7+ -A / 

where [b , b + ] = 1 , [7 , 7 + ]+ = 1 and 7 2 = (7 + ) 2 = . It is straightforward to see that 
the corresponding realization of the loop superalgebra will have only two nonzero com- 
mutators. Hence, one can consider a combination of Gaudin and oscillator realizations 



L(X) = L osc (X) + L(X) , 



(3.47) 
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preserving all the properties mentioned above of the corresponding B -operators with 
appropriate changes in hamiltonians and Bethe equations: /?m(a*j) — * Phdi^j) + fJ>j ■ 

Further modifications can be obtained considering quasi-classical limit of the quan- 
tum spin system with non-periodic boundary conditions and corresponding reflection 
equation. 

The expression of the eigenvector of a solvable model in terms of local variables 
parameterized by sites of the chain or by space coordinates is known as coordinate Bethe 
Ansatz || . The coordinate representation of the Bethe vectors gives explicitly analytical 
dependence on the parameters {^i}^ 1 and {z a }^ useful in a relation to the KZ equation 
(Section 4). Using the Gaudin realization (|3.6|) of the generators 



N + N Y 



a=l a=l 



and the definition of the creation operators ( |3.28| ), one can get the coordinate represen- 
tation of the B -operators: 



N 



S M ( W ,/x 2 , = X>+ • • ■ u+J ff J] Viiv®}?' 1 ,**) , (3-48) 



a=l 



where the first sum is taken over ordered partitions it of the set (1,2, ...,M) into 
subsets K a , a = 1,2, . . . , N , including empty subsets with the constraints 

U/C a = (l,2,...,M) , /C a f|X: 6 = fora/6. 

a 

The corresponding subset of quasimomenta 

(a) _ (a) _ (a) _ e if 

where | /C a | is the cardinality of the subset K, a , and jk < jk+i , entering into the 
coordinate wave function 

f({^m}f ] ;z) = (-l) p(ff) (K(i) -^(2))K(2) -^(3))---h/C| 



K\ 



Due to the alternative sum over permutations a € «Si;n this functions is antisymmetric 
with respect to the quasi-momenta. Finally the first factor in ( 3.48| ) (v^ • • • v^ M ) means 
that for j m G K, a corresponding indices of u+ are equal to a so that i>+ = v+ . One 

can collect these operators into product n^=i W)'^"' , consequently we have an extra 
sign factor (— 1)p( 7F ) . 

This coordinate representation is similar to the representations obtained in 0, ^, ||| 
for the Gaudin models related to the simple Lie algebras. The Zi -grading results in 
extra signs, while the complicated structure of the Bm -operators is connected with the 
fact that (Vj~) 2 = ^ while for j ^ k and anticommute. 
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4 Solutions to the Knizhnik-Zamolodchikov equation 



Correlation functions ip(zi, . . . , z n ) of a two dimensional conformal field theory satisfy 
the Knizhnik-Zamolodchikov equation pO] 



Zb 



(4.1) 



where Y" are generators of an orthonormal basis of a simple Lie algebra in a finite 
dimensional irreducible representation V a and i/j(zx, . . . , z n ) is a function of N complex 

N 

variables taking values in a tensor product ® V a . The operator on the right hand side 

a=l 



of ( f4,ip is a Gaudin hamiltonian (1.1). 

A relation between the Bethe vectors of the Gaudin model related to simple Lie 
algebras and the solutions to the KZ equation is well known for sometime |7|, || . Approach 
used here to obtain solutions to the KZ equation corresponding to Lie superalgebra 
osp(\\2) starting from B- vectors ( 3.3C| ) is based on f|]. 

A solution in question is represented as a contour integral over the variables m - ■ ■ 



(4.2) 



ifj(zi, . . . z N ) = j> . . . j> 4>{fl\z)^{fl\z) dm . . .dfi M , 
where an integrating factor (j){fl\z) is a scalar function 

M N 

= I I (Mi - Mj)" I I ( z * - z b )-^r I fj TT(m* - z c )^ 



M N / M N 

%<j a<b \k=l c=l 



(4.3) 



and i£(j2\z) is a Bethe vector ( 0.3O[) where the Bethe equations are not imposed. 
Differentiating the integrand of (|4.2|) one finally gets 



M 



3=1 



-IV 



(4.4) 



A closed contour integration with respect to /Ui, . . . , \im will cancel the contribution from 
the terms under the sum in ( |4.4j ) and therefore ifi(zi, . . . zn) given by ( |4.2| ) satisfies the 



KZ equation. This follows from the form of the integrating factor (4.3) 



k8 z 



nd. 




E 



M 9 ! 



<t> = (3 M {m)9 , 



the differential identity ( 3.26| ) and commutation relations ( p. 25 ). 

Conjugated Bethe vectors (Bj^fl_)* are entering into solution ip(zi, 
dual KZ equation 



k—^(zi,... ,z N ) 



^( Zl ,...,z N )H^ . 



(4.5) 
(4.6) 

zn) of the 
(4.7) 
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The scalar product yip(zi, . . . , zn) , ■ • • > z n)J does not depend on {zj}^ and its 

quasi-classical limit k — > gives the norm of the Bethe vectors Q due to the fact that the 
stationary points of the contour integrals for k — > are solutions to the Bethe equations 
dS/d/jLj = f3 M (Hj)(p, 

N M N M 

S(fl\z) = K ln0 = ^/ a / b ln(z a - z b ) + ^ lnpix, - ^) - ^^Z a ln(2i a - ^) . (4.8) 

a<fe i<j a=l j=l 

According to the observation in the end of Section 3 analytical properties of the Bethe 
vectors of the osp(l|2) Gaudin model coincide with the analytical properties of the sl(2) 
Gaudin model. Thus, the expression for the norm of the Bethe vectors fl3.3C ) obtained 
as the first term in the asymptotic expansion k — > coincides also 

||^ 1 ,.., M ; {M f)f=det(^|-), (4.9) 
= y _L y 1 d*S _ 1 for 

(4.10) 

Finally we notice that the modification of the Gaudin hamiltonians we mentioned at the 
end of the previous Section can be transfered to the corresponding modification of the 
KZ equations. Both modifications are related with extra factors in the integrating scalar 
function 



5 Conclusion 

The Gaudin model corresponding to the simplest non-trivial Lie superalgebra osp(l\2) 
was studied. A striking similarity between some of the most fundamental characteristics 
of this system and the sl(2) GM was found. Although explicitly constructed creation 
operators Bm (|3.28j ) of the Bethe vectors are complicated polynomials of the generators 
v + {\) and X + {\) , the coordinate form of the eigenfunctions differs only in signs from 
the corresponding states in the case of sl(2) model. Moreover, the eigenvalues and the 
Bethe equations coincide, provided that the sl(2) Gaudin model with integer spins is 
considered. 

Let us point out that using the method proposed in this paper one can construct ex- 
plicitly creation operators of the Gaudin models related to trigonometric Izergin-Korepin 
r-matrix [16] and trigonometric osp(l\2) r-matrix 21]. Similarly to the simple Lie 



algebra case solutions to the KZ equation were constructed from the Bethe vectors using 
algebraic properties of the creation operators Bm and the Gaudin realization of the 
loop superalgebra C+(osp(l\2)) . This interplay enabled us to determine the norm of 



eigenfunctions (4.9) 



The difficult problem of correlation function calculation for general Bethe vectors 

c (Mf ; {Mf; ) = ( «-> b mK • • • m) H h(\ k )B M (^ . . . 



k=l 



was solved nicely for the sl(2) Gaudin model in || using the Gauss factorization of the 
loop algebra group element and the Riemann-Hilbert problem. The study of this problem 
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for the osp(l\2) GM is in progress and the following expression for the scalar product of 
the Bethe states is conjectured (cf. ||) 

(n_ , B* M (v u . . . v m )B m (m,. ..hm)^-) = (-1) P(,T) detM a , 

o-G<S M 

where the sum is over symmetric group S\j arid M X M matrix M. a is given by 



Pi ~ V °V) k^j ^ ~ ~ u °( k )) 

1 

(Pi -Vk)(v a (j) ~ v <r(k)) 



M% = : . lor 7 ./, -J. 2 U. 
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